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' The effective quark Lagrangian is derived from the QCD action by 

. means of averaging over vacuum gluonic fields. In the hmit of large 

\ Nc and keeping for simplicity the lowest order gluonic correlator one 

obtains nonlinear equations for the quark propagator, exhibiting scalar 
confinement. Connection to the quark zero mode density is discussed. 



Abstract 



D , 

• • 1 Introduction 

> ■ 

^ ' The chiral symmetry breaking (CSB) in QCD together with confinement are 

■ two most important nonperturbative phenomena, which define the structure 

of vacuum and hadrons. The simplest situation where both these phenom- 
ena define the dynamics is the system of one light quark and a very heavy 
antiquark, to be considered static source. 

In the case of QED the relativistic system of a light and heavy charges 
was considered among other things in [1] . 

The Feynman-Schwinger representation [2] studied in [1] allows to make 
the infinite mass limit and obtain the effective Dirac equation for the light 
charged particle (this is in contrast to Bethe-Salpeter equation, where the 
one-body limit is not recovered at least in lowest order approximation for 
the kernel). 
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The attempt to extend these resuhs to QCD [3] meets with serious dif- 
ficulties, and one needs the real understanding of dynamics of confinement 
and CSB in their interconnection. 

While confinement was explained recently as being produced by specific 
vacuum correlators of gluonic fields [4] (of the form of correlators of effective 
magnetic monopole currents), the phenomenon of CSB is still not understood. 

In this paper we suggest a systematic method to treat CSB in connection 
with confinement. We start with the QCD Lagrangian and derive from that 
the effective Lagrangian of light quarks, assuming that certain gluon field 
correlators are nonzero, which are known to yield confinement (i.e. linear 
potential) for static quarks. It is not clear from the beginning what will 
happen for light quarks ( with vanishing mass) and whether they would be 
confined at all. 

Our main concern in what follows is to keep Lorentz and gauge invari- 
ance. The effective quark interaction is nonlocal and to have gauge-invariant 
equation one should consider qq Green's functions. 

As mentioned above the simplest setting for which confinement and CSB 
can be studied in the gauge-invariant way, is the problem of a light quark 
propagation in the field of the static antiquark. To simplify matter we start 
with the Gaussian correlators for gluon fields and derive the selfconsistent 
equations for the light quark propagator (with a string effectively connecting 
it to the static source). We show that CSB occurs due to the string (linear 
confinement), which shows up in the fact that effective interaction becomes 
Lorentz scalar. As another evidence of CSB in [5] the chiral condensate is 
computed and shown to be nonzero. The strongly modified and extended 
version of this study in [5] contains also another derivation of scalar confine- 
ment. 

Those results were derived actually for the case of one light quark (quench- 
ing approximation or iVg — > cxd ). To treat the case of two and more light 
quarks the corresponding term in the effective Lagrangian may be studied, 
and one is naturally led to the equation for the Green's function of two 
(or more) light quarks. The study of two and three quark systems can be 
performed in the framework of the same formalism and is now in progress. 

The paper is organized as follows. In the second section the general form 
of the effective Lagrangian for light quarks is given, and the equation for the 
Green's function of one fight quark in the field of the static source is derived. 
In the third section this equation is studied and solutions are derived, which 
exhibit scalar confinement and hence CSB. 
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In the fourth section the connection of the chiral condensate to quark 
zero modes is illustrated in the framework of the present approach. 

The conclusion is devoted to a summary of results and the discussion of 
prospectives. 

2 Derivation of the effective Lagrangian for 
the hght quark 

To make discussion of confinement and chiral symmetry breaking fully gauge 
invariant, we consider the gauge invariant physical amplitude - the Green's 
function 5" of a light quark of mass m in the field of a static antiquark placed 
at the origin. The propagator for the latter in positive time direction can be 
taken as (we consider below Euclidean space-time). 

Sq{T) = exp igJ^A,{f^ 0,r)dr (1) 

The Green's function S^q can be written as an integral 

(2) 

Since SqQ is gauge invariant, one can choose any convenient gauge for 
and our choice will be the modified Fock-Schwinger gauge, namely for any 

point X one introduces the curve C{x,s),s < 1 with the tangent vector t^{s) 
and coordinate z^{x, s). The gauge condition (for some general class of the 
curves C {x, s) considered below) is 

A^,{z)t^ = 0. (3) 

We choose the curve C to go from the point x along some path (to be 
specified later) to the 4-axis, and then run along it, so that A4('r — 0, r) = 0. 
In this gauge Gq reduces to the factor and one can now consider the 
integration over DA in (2) as a statistical averaging process and use the 
cluster expansion [6], {a,b - color indices) 

J £)^e/'^^(^^^(^)'^(^)'^'^-/^^^'.'^*^ exp J ilj^{x)Aab{x)Mx)d^x > 

= exp{ / d^x'4>'^{x)'-^fj, < An> ip{x)+ (4) 
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We have denoted the higher order cumulant contribution in (4) as ... 
and shall disregard it here, for the discussion of these terms the reader is 
referred to [5] . The first term < > vanishes due to the gauge and Lorentz 
invariance of the vacuum, while the second can be expressed through the 
field strength correlators < Fni,{u)Fx^{u') > in the gauge (3) as follows 



\ ( \ f ^ dZa[s,x , .dzp 



The representation (5) satisfies condition (3) for the class of curves C 
considered below. Using (5) one can rewrite the average < AA > in (4) as 

^ / V /V du0 dvs 
X < K,(u)F,,(v) > (6) 

where ab, cd are fundamental color indices. The gauge-nonivariance of the 
correlator < FF > in (6) is only apparent and one can introduce a factor, 
equal to unity in the gauge (3), which makes the correlator explicitly gauge- 
invariant, namely 

< F{u)F{v) >=< F{u)-^F{v)-^+ > (7) 
where ^ is the product of 2 parallel transporters 

* = ^ciuMi^) (8) 

and 

^c{u) = P exp ig dz^A^ (9) 

where the path-ordered contour integral $ is taken along the contour C. 

Since in the gauge (3) one has — = 1, we shall below omit those fac- 
tors. For the correlator < FF > one can use the parametrization suggested 
in [7] 

g'^ < Fi^{u)Fi,^,{u ) >ab= Sab{Sii'S^^> - Sif,'6i'^)D{u - u') + A^^^ (10) 
where A*^^) is proportional to a full derivative, its exact form is given in [7]. 
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In what follows we shall consider only the term D in (10) since it con- 
tributes to the string tension, while A*^^^ does not. Namely using (10) it was 
obtained in [7] that the string tension a - the coefficient in the area law of 
the Wilson loop, < W{C) >= exp{—aarea) is equal to 



1 /"OO 
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d\D{u) (11) 



Thus a characterizes the confinement of static quarks, and our goal is to 
understand how the dynamics of light quarks is expressible through a and 
whether a correlates with CSB. 

Keeping only the term D in (10) and neglecting higher order correlators 
hke < AAA >, one obtains in (4) the following effective Lagrangian for the 
light quark: 

>Ce//('0'^'0) — j 'il^~^{x){—id — im)'ip{x)d'^x+ 

^ I d'xd'y{^p^{x)^,U^mt{y)i,'i^a{y))^ (12) 

^J^l^i'{x,y) 

where we have defined 

dun dvs 



rz rw I 

Jnn'{z,w) = dUa dv^{Sa'YSps " SasSfi-y)- 



■X 



xD{u-v) (13) 

In what follows we disregard the perturbative contributions to C^ff, since 
they have nothing to do with CSB. The mass m is supposed to be defined 
at the typical hadronic scale of 1 GeV, and we shall not be interested in its 
evolution to lower scales. 

Prom the effective Lagrangian (12) in the large Nc hmit one can easily 
derive the equation of Dyson-Schwinger type for the selfenergy part, which 
we shall denote by M and the qQ Green's function S. (We denote by S 
the form (2) where instead of Sq only the parallel transporter is retained 
(cf. Eq.(l))). The definition and the diagrammatic representation of M is 
done in the same way, as in the NJL model [8], since the structure of the 
Lagrangian (12) is similar to that of NJL however nonlocal. 

The main essential difference is the presence of the string , connecting 
the light quark to the static source, this part is concealed in J, Eq.(13) and 
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therefore the selfenergy part M is actually not the set of the one-particle- 
irreducible diagrams, but rather the qQ interaction kernel. 

In the configuration space the equations for M and S are readily obtained 
from (12) noting that in the large Nc limit one has to replace a pair of ip, 
operators in (12) as 

Mx)i^b{y) ^< Mx)^biy) NcS{x,y), (u) 

and finally one obtains 

iM{z,w) = Ji^„{z,w)'yf,S{z,w)'y^ (15) 

{-id, - im)S{z, w)-i j M{z, z')S{z', w)d^z' = 5^^^ {z - w) (16) 

The system of equations (15-16) is exact in the large limit, when higher 
correlators are neglected and defines unambiguously both the interaction 
kernel M and the Green's function S. One should stress at this point again 
that both S and M are not the one-particle operators but rather two-particle 
operators, with the role of the second particle played by the static source. It 
is due to this property, that S and M are gauge invariant operators, which 
is very important to take confinement into account properly. Had we worked 
with one-particle operators, as is the habit in QED and sometimes also in 
QCD, then we would immediately loose the gauge invariance and the string, 
and hence confinement. 

3 Chiral symmetry breaking by the solutions 
of equations (15), (16) 

As was discussed in the Introduction, CSB may reveal itself in the heavy-light 
system in two ways: i) as the presence of the scalar part in the effective mass 
operator M and ii) as the nonzero quark condensate tr5'(0, 0) ~< >. 

Both points have been studied in [5] using the relativistic WKB method 
and the presence of scalar confinement and nonzero quark condensate was 
shown as a consequence of the selfconsistent WKB solution. 

In this paper we shall give another argument in favour of the scalar con- 
finement, using an approximate asymptotic form of the solution of (15-16) 
at large distances. 
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The idea which enables us to find this approximate solution is based on 
the following simple physical picture. At large distances the string formed 
between the quark and the static source becomes heavy due to the length of 
the string even for the massless quark. Therefore the situation is similar to 
the case of the heavy quark, where the exact solution is known in the limit of 
the infinitely large mass, and one can use an expansion in powers of inverse 
mass of the string. 

To introduce this expansion we first specify the kernel J^^ in (15) and to 
this end define the contour C in (5) and the form of the correlator D[x) in 
(10). 

For the latter it is convenient to choose the Gaussian form, 

D{x) ^ D,exp{-^1 ^'^2^2 (17) 

where the use of equation (11) was made. 

To calculate J^ui^, w) we choose the special case of the gauge (5), which 
may be called "the planar gauge": the contour C{x) is formed by the three 
straight-line pieces: i) first a perpendicular to some "gauge" plane, passing 
through the 4-th axis, then ii) a perpendicular to the 4-th axis, and finally 
iii) a cut along the 4-th axis to 0:4 = —00. One can choose the "gauge plane" 
passing in the middle between the points z and w in J^y{z4^ — W4, z,w). To 
simplify calculations let this plane be the 14 plane, and then omitting for the 
moment the contribution of the cut i) in the integral (5) (we shall discuss 
its contribution below) one obtains the following representation of the only 
nonzero component of J/^i, = J44 

2 11 2 

_ (24-1114) . _(zis-wit) 

Mz,w) = e ~^a^ J J dsdte (18) 

f 

In the limiting case when 3- vectors z and w coincide and lie in the gauge 
plane, the expression (18) simplifies 

MzM - ^expi-^-^^^) (19) 
V^Tg 4T/ 

We now turn to equations (15 - 16) and look for the solution in the form 
of an expansion 

S = 50 + ^1 + ... 
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(20) 

M = M0 + M1 + ... 

with 

Mo = J44745'o74 (21) 

and 

h 

So(z,w) = '-{e{h){l + ^4)exp[- J M{h',z,z)dh']^i^\z,w)+ 



h 

+ e{-h){l - ^^)exp[J M{h',z,z)dh']ipl^\z,w)} (22) 


Here notations are used 

h = Z4-W4, rf^ = (fi^o^ = 5(3) {z - w) (23) 

and M is to be defined from tlie solution of (15 - 16). 
Acting with the operator — i J^74 on (22) one has 

- i(9474^o(/i, z, w) = 6^'^\z -w)+ iM(/i4, z, z)So{h, z, w) (24) 
From (21) one obtains 

Mo{h,z,w)^ ^^d^^\z-w)e~^N (25) 

where 

h h 

N = e{h){l + 74)ea;p(- J Mdh') + e{-K){l - ^A)exp{ J Mdh') (26) 



To simplify matter further wc take the so-called string limit, i.e. tend 
Tg to zero while keeping string tension a constant. In this case it will be 
sufficient to consider M independent of time h! in (22), and one obtains 

j d^z'Moiz, z')So{z', w) = a\z\So{z, w) (27) 
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In derivation of (27) the condition cr|5|Tg <^ 1 was used, i.e. the string 
hmit is taken before any other asymptotics is considered. 

The space derivatives in (16) can be shown to give subleading terms, as 
compared to (24) and (27), this is exactly the same situation as in the heavy 
mass limit considered in [5]. Indeed, inserting (24) and (27) into (16) on has 

M{z,z) = a\z\+m (28) 

— * 

At the same time the spacial derivative ^dSo does not contain the large 
parameter — > oo, and therefore should be taken into account in the next 
order of the expansion (20). The appearence of the derivative ^d5^^\z — w) 
signals that the local approximation (23) used for 5*0 should be in next orders 
replaced by the nonlocal form, which was obtained in [5], 

Mr^S{x,y) ~ -,^^K^{a^\x-y\)5{l-cose^y) (29) 

where Kq is the McDonald function. The r.h.s. of (29) is the smeared 
5-function at large \x\ ~ which should replace (23) in the next approxi- 
mation. 

One might improve the initial approximation (22) keeping in (16) the local 
form of M (25), (since it enters under the integral) but replacing ^Pq\(^q'^ 
in (22) by some nonlocal form. Instead of writing explicit formulas, one can 
obtain the simple static form of equation (16) if one integrates both sides of 
(16) over all infinite interval oi h — z^ — w^. One obtaines 

[-i^d - i(m + a\z\)]S{z,w) = S^^\z - w) (30) 

where we have defined 

/oo 
dhS{h, z, w) (31) 
-oo 

Both equations (28) and (30) signify the scalar confinement at large dis- 
tances, independently of the value of the quark mass m. 

This statement agrees with the conclusion of the paper [5], where another, 
less direct method was used. 

At this point one can compare our results with those in [9] and [10], 
where the limit of the heavy quark mass m was used (the authors of [9] have 
exploited the form (15) with S replaced by 5*0 as the basic equation), and 
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at the same time to elucidate the dependence of results on the choice of the 

gauge. 

One should stress first that the gauges used in [5], [9] and [10] on one side 
and in the present paper are different. 

This fact does not play any role if one considers the total effective action, 
i.e. the form (12) with addition of all higher order correlators as in [5], Eq. 
(163). However the Gaussian approximation for L^/f and M, as in (15) 
depends on the gauge, or (if one considers the quark trajectory L with the 
line C{x) attached at each point of L) on the shape of the gauge surface, 
this dependence being cancelled by higher order correlators. To minimize 
the influence of higher-order correlators one should choose the shape of the 
gauge surface as close to the real physical world sheet of the string as possible, 
and the latter is known to be the minimal surface. 

The situation here is parallel to that of the perturbative QCD series, 
where keeping a few first terms one obtains a dependence on the normal- 
ization mass II, which should not enter into the physical result. The usual 
prescription is to choose /i in such a way as to minimize the neglected terms 
of the perturbative series and the recipe is to take /j, of the order of the inverse 
size of the system. 

In accordance with what was said above it is profitable to choose the gauge 
surface of minimal area, and our choice of the plane instead of the surface in 
[5], [9] and [10] strongly reduces the contribution of the spacial projections of 
the gauge surface. The latter enter in the final result being multiplied by the 
color-magnetic correlators < Fij{z)Fim{w) >, i,j,l,m = 1,2,3. Comparison 
of (28) with the results of [5], [9] and [10] shows that (28) contains only 
the contribution of the color-electroc correlator which agrees with that of 
[5], [9] and [10], while the latter papers contain also the contribution of 
magnetic correlators. Hence the color-electric contribution is surface (and 
gauge-) independent, while the color-magnetic one depends on the shape of 
the gauge and can be excluded for the best (mininal) choice of the gauge 
surface (the above is true only for the zero orbital momentum, otherwise the 
minimal surface is close to the helycoid, rather than to the plane, and the 
magnetic contribution survives). 
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4 Quark zero modes vs field correlators 



We have seen in previous sections that the nonhnear equations (15), (16) give 
rise to the phenomenon of CSB, which reveals itself in our problem in two 
ways: i) it provides scalar confining interaction for the light quark ii) there 
appears a standard chiral condensate < i/jiIj >. 

A natural question arises at this point: a folklore understanding of CSB 
is that it is due to quasizero global quark modes in the vacuum. An exact 
relation [11] exists, which connects chiral condensate to the density z/(A) of 
quasizero modes in the vacuum at A < m, and in the chiral limit (m — > 0) 
one has [11] 

< >= (32) 

Here A is an eigenvalue of the 4d Euclidean equation for the quark in the 
vacuum field 

- iD'il)n{x) = \ni^n{x) (33) 

The density v{\)d\ is the averaged over all fields {A^} number of the states 
A„ per interval d\. 

It is a popular belief that the quasizero modes necessary for CSB due 
to (32) are descendant from the local zero modes on the topological charges 
(instantons or dyons), and their density is therefore proportional to the den- 
sity of instantons (dyons). There are the instanton model [12] and the dyon 
model [13] of the QCD vacuum, which explain CSB in this way. 

Whether these models are realistic or not, is the open question, but the 
Banks-Casher relation (32) holds independently of that, and if the method of 
the present paper proves CSB due to the field correlators (even in the Gaus- 
sian approximation) , one should explain the origin of the quasizero modes in 
(32). 

To do this we consider first the case of Abelian fields. As was stressed 
above in the paper, CSB is due to the correlator D{x), and the latter in 
the Abelian case can be connected to the correlator of magnetic monopole 
currents [7,4] 

< hi^&w >= (~jf. - - y) (34) 

In the nonabclian case one can use the Abelian projection method (APM) 
[14] , to separate out of the field ^4^ and field strength F^,^ the monopole and 
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photon part, and the part of "charged gluons" . The latter contributes around 
10% to the effective action. In this case one can connect the monopole current 
obtained by APM with the correlator D as in (34). 

Now each magnetic monopole with the world line of the length T produces 
T/h quasizero modes, where h is the size of the monopole [15]. The interval of 
modes produced by monopoles is AA, threfore the 3d density of monopoles 

multiplied with 3d volume V3 and with T/h gives the total number of 
quasizero modes in the intervals AA, i.e. z/AA = ^V^n^- 

Taking into account that V4 = TV3, one has a relation ~ ^^(AA6). 
The last factor (AA6) is equal to unity if the quasizero modes of monopoles 
are mixed due to interaction into the whole interval A A ~ 1/6. (As will 
be seen under this assumption one obtains the estimate (35) for the chiral 
condensate.) On the other hand one can estimate the 3d density of magnetic 
monopoles, from < ji3{x,x^)jg{{)^Xi) > integrating over d^x. (The correlator 
< j{x)j{0) > estimates probability of finding a monopole at the point x, 
if there is one ai x — 0. Integrating over d^x one finds the probability of 
having a monopole at x — 0, while another is anywhere; fixing X4 means that 
the probability refers to a given moment. Note the condition on magnetic 
charge / j^^d^x = 1, which yields the correct normalization.) Hence one gets 
an order of magnitude relation 



where we have assumed for D{x) the form D{x) — D{0)f{^) and f{y) is an 
exponential or Gaussian with /(O) = 1. 

Finally, taking into account that -D(O) ^ g'^ < F(0)F(0) > one obtains 



This estimate coincides with our result obtained from the quasiclassical cal- 
culation in [5]. Numerically (36) is - (SOOMeV)^, i.e. a reasonable order 
of magnitude. Thus the very existence of the "wrong" correlator D{x), vio- 
lating the Abelian Bianchi identity may bring about monopole currents and 
ascociated with those zero modes. 




(35) 



< Ipljj >Ri 



< F^{Q) > Tg 



(36) 
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5 Conclusion 



The main question posed at the beginning of the paper -what is the exphcit 
form of dynamics, i.e. equations for a hght quark in the field of a heavy static 
source. The answer to this question is given by equations (15), (16). 

Those are Dirac-type equations with nonlocal nonlinear interaction. We 
have shown that these equations allow for solutions which yield scalar con- 
fining interaction. Therefore at large (as compared to Tg) distances one deals 
effectively with a local Dirac equation with a scalar potential. This is in con- 
trast to the QED case considered in [1] , where the one-body limit yields the 
well known Dirac equation with vector Coulomb potential. The local scalar 
potential was assumed for the heavy-light qq system in [3] , where numerical 
and analytic results for spectra and wave functions are presented. 

The appearance of the scalar interaction provides (together with nonzero 
chiral condensate) an evidence for CSB, and the fact that these effects are 
proportional to the string tension, tells about interconnection of CSB and 
confinement. There are several direction for the extension of this study. 
First, one should consider the case of 2 and 3 light quarks and derive the cor- 
responding chiral Lagrangian. Second, with CSB at hand one can attack the 
problem of constituent quark mass. These topics are planned for subsequent 
publications. 

A large part of results of this paper have been obtained when the au- 
thor was a guest of the Institute of Theoretical Physics in Utrecht. The 
cordial hospitality and interesting discussions with Professor John Tjon are 
gratefully acknowledged. 
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